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Abstract 

We employ monomer-resolved Molecular Dynamics simulations and theoretical considerations to 
analyze the conformations of multiarm polyelectrolyte stars close to planar, uncharged walls. We 
identify three mechanisms that contribute to the emergence of a repulsive star-wall force, namely: 
the confinement of the counterions that are trapped in the star interior, the increase in electrostatic 
energy due to confinement as well as a novel mechanism arising from the compression of the stiff 
polyelectrolyte rods approaching the wall. The latter is not present in the case of interaction 
between two polyelectrolyte stars and is a direct consequence of the impenetrable character of the 
planar wall. 

PACS numbers: 82.70.-y, 82.35.Rs, 61.20.-p 
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I. INTRODUCTION 



Polyelectrolytes (PE's) are polymer chains that carry ionizable groups along their back- 
bones. Upon solution in aqueous solvents, these groups dissociate, leaving behind a charged 
chain and the corresponding counterions in the solution. Chemically anchoring / such PE- 
chains on a colloidal particle of radius -Rd gives rise to a spherical polyelectrolyte brush 
(SPB). In the limit in which the height R of the brush greatly exceeds R^, one talks in- 
stead about polyelectrolyte stars. A great deal of theoretica l^i^i^i'^i^i^i'^i^i^i^^i^'^i^^i^^i'^^i'^^i'^^i'^''' 
and experimenta l^^i^^i^^i^^i^^i^^i^^i^^i^^i^''' effort has been devoted in the recent past with the 
goal of understanding the conformations and the interactions of SPB's and PE-stars. The 
reasons are many-fold. The grafted PE-chains can provide an electrosteric barrier against 
flocculation of the colloidal particles on which the chains are grafted, rendering the systems 
very interesting from the point of view of colloidal stabilization. Moreover, they can act as 
control agents for gelation, lubrication and flow behavior. From the point of view of funda- 
mental research, PE-stars are a novel system of colloidal particles that combine aspects from 
two different parts of colloid science: polymer theory and the theory of charged suspensions. 
Accordingly, the derivation of the effective interaction between PE-stars has led to their 
description as ultrasoft colloids and to theoretical predictions on their structural and phase 
behavior-^^ that have received already partial experimental confirmation.^^ Moreover, cer- 
tain similarities in the structural and phase behavior of PE-stars and ionic microgels have 
been established, demonstrating the close relationship between the two systems.— >2S 

By and large, the theoretical investigations involving PE-stars have been limited to the 
study of either single PE-stars or bulk solutions of the same. However, a new field of 
promising applications and intriguing physics is arising when PE-stars or SPB's are mixed 
with hard colloids or brought in contact with planar walls. Indeed, PE-stars can be used 
to model cell adhesion and are also efficient drug-delivery and protein-encapsulation and 
immobilization agents . On the other hand, hydrogels, which are physically similar to 
PE-stars, adsorbed on planar walls, form arrays of dynamically tunable, photoswitchable or 
bioresponsive microlenses.— i2S»2L2& The encapsulation properties as well as the characteristics 
of the microlenses depend sensitively on the interactions between the PE-stars and the 
colloidal particles or the wall, respectively. Therefore, there exists a need to undertake a 
systematic effort in trying to understand these interactions physically and make quantitative 
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predictions about the ways to influence them externally. In this paper, we take a first 
step in this direction by considering a PE-star in the neighborhood of a planar, uncharged 
and purely repulsive wall. We analyze the mechanisms that give rise to an effective star- 
wall repulsion and identify the counterion entropy and the chain compression against the 
planar wall as the major factors contributing to this force. Our findings are corroborated 
by comparisons with Molecular Dynamics simulation results for a wide range of parameters 
that structurally characterize the stars. This work provides the foundation for examining 
effects of wall curvature and charge. 

Our paper is organized as follows: Sec. HHis devoted to the description of the system we 
are investigating, the simulation model, and the simulation techniques used. Moreover, we 
specify the physical quantities of interest. Our theoretical approach is addressed in detail in 
Sec. mil In Sec. IIVI we quantitatively compare and discuss the respective results. Finally, 
we conclude and present an outlook to possible future work in Sec. El 

II. SIMULATION MODEL 

We start with a definition of the system under consideration including its relevant param- 
eters and a description of the simulation model used. We study a dilute, salt-free solution of 
PE-stars confined between two hard walls parallel to the x-y-pla.ne at positions z = ±t/2, 
resulting in an overall wall-to-wall separation r. We assume a good solvent that is only 
implicitly taken into account via its relative dielectric permittivity e ~ 80, i.e., we are deal- 
ing with an aqueous solution. To avoid the appearance of image charges , we assume 
the dielectric constants to be the same on both sides of the respective confining walls. It 
turns out, however, that this is not a severe assumption, because the effective charge of a 
PE-star is drastically reduced compared to its bare value, due to the strong absorption of 
neutralizing counterions.-^'^'S^ Hence, the influence of image charges can be expected to be 
of minor importance. 

The PE-stars themselves consist of / PE-chains, all attached to a common colloidal core 
of radius R^, whose size equals the monomer size and is therefore much smaller than the 
typical center-to-end length of the chains for all parameter combinations. The introduction 
of such a core particle is necessary to place the arms in the vicinity of the center, where 
the monomer density can take very high values. The theoretical approach pertains to the 
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limit of vanishingly small core size. In order to remove effects arising from the small (but 
finite) value of the core in the simulation model and to provide a comparison with theory, we 
will henceforth employ a consistent small shift of the simulation data by R^- This standard 
approach has been applied already, e.g., to the measurement of the interactions between 
neutral^ and charged stars^^, in order to isolate the direct core-core interaction effects from 
the scaling laws valid at small interstar separations. 

The PE-chains are modeled as bead-spring chains of Lennard- Jones (LJ) particles. 
This approach was first used in investigations of neutral polymer chains and stars^^*^^^ and 
turned out to be reasonable. The method was also already successfully applied in the case of 
polymer-colloid mixtures^^ or polyelectrolyte systems.— i^*^ To mimic the above mentioned 
good solvent conditions, a shifted and truncated LJ potential is introduced to depict the 
purely repulsive excluded volume interaction between the monomers: 



^Lj(r) 



Ae 



LJ 



(r)"-(^r + i 



(1) 

r > 



Here, r is the spatial distance of two interacting particles. The quantities ulj and ^lj set 
the basic length and energy scales for the system, respectively. In what follows, we fix the 
system's temperature to the value T = 1.2ei^j/k-B, where denotes Boltzmann's constant. 
The PE-chain connectivity is modeled by employing a standard finite extension nonlinear 
elastic (FENE) potential>^^ 

f fepENE f Rq \ ]y. 



r<Ro 

(2) 



oo r > i?o, 

with a spring constant fcpENE = 7.0£lj- The divergence length Rq limits the maximum 
relative displacement of two neighboring monomers and is set to Rq = 2.0crLj in the scope 
of this work. The described set of parameters determines the equilibrium bond length, in 
our case resulting in a value Iq = 0.97(Tlj. 

When modeling the interactions between the monomers and a star's colloidal core, the 
finite radius -Rd of the latter has to be taken into account. All monomers experience a 
repulsive interaction with the central particle, in analogy to Eq. reading as 

oo r < Rd 

v^j = <; (3) 

^Lj(r-i?d) r>i?d. 



In addition, there is an attraction between the innermost monomers in the arms' chain 
sequence and the core which is of FENE type and can be written as (cf. Eq. (j2I)): 

oo r < Rd 

V^ENE = { (4) 

Vfene('^ - ^d) r> Rd- 

The chains are charged in a periodic fashion by a fraction a in such a way that every 
(l/Q;)-th monomer carries a monovalent charge. Consequently, there is a total number of 
Nc = otfN monomer ions per PE-star. To ensure electroneutrality of the system as a 
whole, we include the same amount of oppositely charged, monovalent counterions in our 
considerations. Since the latter are able to freely move, they have to be simulated explicitly. 
Furthermore, they are of particular importance because they are expected to crucially affect 
the physics of the system. One example is the aforementioned fact that they induce a 
reduction of the stars' bare total charges. 

Two charged beads with spatial distance r^j interact by a full Coulomb potential, i.e., the 
electrostatic interaction energy is 

PVcon^in,) = = Ab^. (5) 

Thereby, = ±1 are the valencies of monomer ions and counterions, respectively, e 

denotes the absolute value of the unit charge, and [3 = l/k-oT is the inverse temperature. 
In the above equation, the so-called Bjerrum length 

Ab = — (6) 

was introduced. It is defined as the distance at which the electrostatic energy equals the 
thermal energy. Thus, it characterizes the interaction strength of the Coulomb coupling. In 
the case of water at room temperature, one obtains Ab = 7. 1 A. In our simulations, we fix 
the Bjerrum length to A = S.Octlj, thus corresponding to an experimental particle diameter 
a = 2.4 A. This is a realistic value for typical polyelectrolytes.— 

For the purpose of completing the set of interaction potentials needed to describe the 
system at hand, we have to define particle-wall interactions. For technical reasons (see 
below), we do not regard the walls as true hard walls, as common in Monte-Carlo (MC) 
simulation studies.— Following the course of our above modeling and based on Eq. ((T)), 
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we assume them to be of truncated-and-shifted-LJ type instead, leading to the following 
monomer-wall interaction: 

^u(^) = ^Ljg-^)+^Ljg+^), (7) 

whereas z refers to the ^-component of the position vector of the particular bead. Likewise, 
the potential function for a star's core interacting with the confining walls is yielded by 
combining Eqs. (jHl) and 

VIT(^) = ^Lj - i?d - ^) + lij - i?d + ^) . (8) 

To finalize this section, we present a short summary of the simulation techniques used. We 
perform monomer-resolved Molecular Dynamics (MD) simulations in the canonical ensemble, 
employing a rectangular simulation box of total volume V = M'^t that contains a single PE- 
star. We apply periodic boundary conditions in the x- and y-directions, while the box 
is confined with respect to the z-direction. Here, we always fix M = r = 120(Tlj. This 
provides a sufficiently large simulation box to suppress any undesirable side-effects a priori 
and emulates a dilute solution of PE-stars (cf. above description of the physical problem 
under investigation). For the numerical integration of the equations of motion, we adopt 
a so-called Verlet algorithm in its velocity form .^' - 'i^'^^i^^ In order to stabilize the system's 
temperature, we make use of a Langevin thermostali^^i^i*^ that introduces additional friction 
and random forces with appropriately balanced, temperature-dependent amplitudes. Due 
to the periodic boundary conditions and the long-range character of the Coulombic forces, 
a straight-forward calculation of the latter pursuant to Eq. is not feasible. Therefore, we 
have to evaluate the forces using Lekner's well-established summation method^^ in its version 
for quasi two-dimensional geometries. Thereby, the convergence properties of the sums 
occurring during the computation are enhanced by a mathematically accurate rewriting, 
allowing a proper cut-off. For performance reasons, the forces have to be tabulated. 

The typical time step is At = 0.002to, with to = (^'^lj/^lj)^''^ being the associated time 
unit and m the monomer mass. In our case, the counterions are taken to have the same 
mass and size as the (charged) monomers. 

In our simulations, we measure the effective core-wall forces as a function of the center- 
to-wall distance D. Furthermore, we examine various other static quantities and their in- 
dependence, namely the center-to-end distance R, the radius of gyration Rg, the number 
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of condensed counterions and the density profiles of all particle species involved. For this 
purpose, for every particular value of D the system is equilibrated for about 5 x 10^ time steps. 
After this equilibration period, we perform production runs lasting between 1 x 10^ — 2 x 10^ 
time steps. We carry out the above described measurements for a variety of arm numbers 
/ (/ =10, 18, 30) and charging fractions a {a =1/5, 1/4, 1/3) in order to make systematic 
predictions for the /- and a-dependencies of all theoretical parameters. In all cases, we fix 
the degree of polymerization, i.e., the number of monomers per arm, to a value = 50. 



III. THEORY OF THE EFFECTIVE PE-STAR WALL INTERACTION 

A PE-star in the neighborhood of a planar, impenetrable wall, undergoes conformational 
changes that modify its (free) energy in comparison to the value it has when the wall is 
absent or far away from the star. This separation-dependent difference between the free 
energies is precisely the effective interaction between the PE-star and the wall. There are 
three distinct mechanisms that give rise to an effective interaction in our case: the change 
in the electrostatic energy of the star, the change in the counterion entropy arising from the 
presence of the geometric confinement as well as contributions from compression of the stiff 
PE-chains against the flat wall. The former two are intricately related to each other, as the 
counterion distribution is dictated by the strong Coulomb interactions, hence they will be 
examined jointly, the latter is a distinct phenomenon arising in the presence of impenetrable 
walls. 



A. The electrostatic and entropic contributions 



To obtain theoretical predictions for the electrostatic-entropic contribution to the free 



energy for a PE-star close to a flat, hard wal 
by and akin to that developed in Refs. 



we employ a mean-fleld approach, inspired 



and 
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Let ps = Ns/Q be the density of 



the solution of PE-stars, where Ng denotes the total number of PE-stars in the macroscopic 
volume Q. In general, a single PE-star of total charge Q is envisioned as a spherical object 
of radius R, embedded in a spherical Wigner-Seitz cell of radius i?w Additionally, the latter 
contains Ac counterions, restricted to move in the Wigner-Seitz cell only and forming an 
oppositely charged background of total charge —Q. The cell's radius is connected to the 
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star density via -Rw = (4vrps/3)^/^. In the case of center-to- wall separations smaller than 
i?w and/or R, we simply cut the counterion cloud and/or the star itself at the confining 
wall and treat them as chopped spheres instead of full spheres. Obviously, the compound 
system composed of a PE-star and its counterions is electroneutral as a whole. Since we are 
interested in the case of dilute P13-solutions only, we can limit ourselves to the consideration 
of a single PE-star. Clearly, for center-to-wall distances D > i?w there is no interaction 
between a PE-star and the wall within the framework of our theory. Accordingly, we will 
deal with the case D < i?w only. Fig. [T] sketches the physical situation and visualizes the 
decisive length scales of the problem. 

A mechanism crucially influencing the physics of PE-systems is the so-called Manning 
condensation of counterions . ^^i^^i^^i^'^'i^^i^^i^^i^'^ A parameter indicating whether or not such 
a condensation effect will occur is the dimensionless fraction ^ = X^Na/R. If it exceeds 
unity, counterions will condense on the arms of the stars.-^^*^^^ In our case, this condition is 
true for all parameter combinations examined (see below). Hence, the condensation effect 
has to be taken into account in our theoretical modeling. For that reason, we partition 
the counteno,. ,„ t.n.ee difeent state, an ansat. a^ead, put fo.,a.d Refs.fe and 
S : Ni of the Nc counterion are in the condensed state, i.e., they are confined in imaginary 
hollow tubes of outer radius Ab and inner radius ctlj around the arms of the star. The total 
volume^ accessible for the counterions in this state is Qi = 7r(A| — a^j)Rf. A number of A''2 
counterions is considered to be trapped in the interior of the star. These ions are allowed to 
explore the overall volume fijn = Qcs{D,R) of the (possibly chopped) star except the tubes 
introduced above, i.e., fl2 = ^cs{D,R) — Qi. Thereby, the volume Qcs of a chopped sphere 
can be derived by straightforward calculations. With Q{R) = AttR^/3, one obtains: 



Note that the total volume of the star is big compared to the volume of the tubes surrounding 
the arms, Qcs{D,R) ^ Qi. Thus, we will make use of the approximation ^2 ~ Qcs{D,R) 
in all steps to come. The remaining = — Ni — N2 counterions are able to freely 
move within the outer shell of volume Qout = ^3 = ^cs(-D, -Rw) — ^cs{D,R) surrounding 
the star, where i7cs(-D,-Rw) is the total volume of the Wigner-Seitz cell. The subdivision 
of counterions in the three different states is also depicted in Fig. ^ showing a star with 5 




(9) 
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arms that are assumed to be fully stretched for demonstration. 

The electrostatic and entropic part Vcc{D) of the effective interaction Ves{D) between the 
PE-star and the wall, kept at center-to-wall distance D, results after taking a canonical trace 
over all but the star center degree of freedom. With J-'cc{D] R, {Ni}) being the variational 
Helmholtz free energy for a system where a single star faces the confining wall, and which 
contains just the electrostatic and entropic parts, it is defined as:— 

VeeiD) = min J^UD; R, {Ni}). (10) 

R,{Ni} 

In principle, the equilibrium values of R and {Ni} are determined by the above minimization. 
We will discuss this point in more detail shortly. Note that we neglected a second, D- 
independent term on the right-hand side of Eq. ()10|) representing the contribution to the 
free energy for an infinitely large center-to-wall separation. Since it makes up a constant 
energy shift only, it does not influence the effective forces between the star and the wall 
we are mainly interested in. The electrostatic-entropic force contributions are obtained by 
differentiating with respect to D: 

Now, we derive concrete expressions for the terms of which J-'ee{D] R, {Ni}) is built up, 
whereby we want to include both the counterions' entropic contributions and the electrostatic 
energy. To keep our theory as simple as possible, we will omit other thinkable contributions, 
like elastic energies^ of the chains or Flory-type terms arising through self-avoidance.-^*^. 
As we will see in Sec. IIV| our theoretical results in combination with chain-compression 
terms to be deduced in what follows, are capable of producing very good agreement between 
the theory and the corresponding simulational data; these simplifications are therefore a 
^ostenoh justified. Consequently, the free energy jFgg(Z); R, {Ni}) reads as: 

3 

J^UD; R, {N,}) = [/es + 5^ S,. (12) 

i=l 

The electrostatic mean-field energy Ucs is assumed to be given by a so-called Hartree-type 
expression. Let pm(i") and Pi{r) denote the number densities of the monomers and the three 
different counterion species, respectively, measured with respect to the star's geometrical 
center. These density profiles then determine the overall local charge density g{r) of our 
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model system. Therewith, we have: 

U^^l-Jji'rSr'fi^. (13) 

It is convenient to separate the total charge density g{r) into two contributions: ftn(r) in the 
interior of the star, i.e., the volume fiin, and foutlr) in the outer region, i.e., the volume flont- 
Let $in(i") and $out(r) be the contribution to the electrostatic potential at an arbitrary point 
r in space caused by the respective charge density. Using this definitions, we can rewrite 
Eq. (O and get 

Ues =l{ [ dV [$i„(r) + <l>o,t(r)] ftn(r) 

+ f dV [<l>i,(r) + <l>out(r)]fout(r)|. (14) 
On purely dimensional grounds, we expect a result having the general form 

l3Ues =—p— ■ h 



R \ R R 



R 

~l~ ^out— out 



"'in— in I t-. I ~r -^"-in— out 



(15) 



R ' R^ 

where we introduced dimensionless functions h^-p arising from the integrations of the prod- 
ucts $Q,(r)^^(r) in Eq. ()14j) . Here, one should remember that is the number of uncom- 
pensated charges of the star and characterizes its effective valency. The specific shape of the 
/i-functions depends on the underlying charge distributions ^in(r) and f?out(r) alone. 

The terms Si represent the ideal entropic free energy contributions of the different coun- 
terion species. They take the form 



(3S.= [ dVp,(r) [\og{p,{T)al)-l] 



+ 3iVaog(^A^, (16) 

where A is the thermal de-Broglie wavelength. In writing the sum in Eq. (|12p . the respective 
last terms in Eq. yield an trivial additive constant only, namely SA^^c log(A/crLj), which 
will be left out in what follows. 

Now, we have to quantitatively address the electrostatic and entropic terms pursuant 
to Eqs. (fT^ to (fT3|) and (fTB|. For that purpose, we first of all need to specify the above 
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introduced number densities,— Pm(r) and Pi{r). Here, we model the arms of the PE-star to 
be fully stretched, or to put it in other words, the monomer density profile inside the star to 
fall of as^*^^ Pm(r) = Pm('") ~ f^ "^- This is a good approximation, as measurements yield 
a scaling behavior with an only somewhat smaller exponent = —1.8, indicating an almost 
complete stretching of the chains.— i^^^^iSiSSi^ Since the monomer ions are placed on the 
backbone of the chains in a periodical manner (cf. Sec. |n}, their density within the interior 
of the star must obviously show an identical r-dependence. Moreover, the profile for the 
trapped counterions exhibits the same scaling, due to the system's tendency to achieve local 
electroneutrality.— Therefore, it seems to be a good choice to use an ansatz p2{r) = A/r"^ 
for the trapped counterions and Qin{T) = Bjr^ for the overall charge density in the inner 
region. Clearly, these distributions have to be normalized with respect to the total number 
of trapped counterions and the effective charge of the star Q*/e = [Nc — Ni — N2) = N3 
by integrating over the related volumes, Q2 and f2in. In doing so, we obtain A = N2/ {2'kRC) 
and B = Q*/{2nRC), with 

(17) 

We presume the condensed counterions to be uniformly distributed within the tubes sur- 
rounding the PE-chains, i.e., we use pi = Ni/Qi. This approach is supported by simulation 
results on single PE-chaina^ and was successfully put forward in previous studies of PE-star 
systems.-^*^ In a similar fashion, we assume an also uniform distribution of the free counte- 
rions within the outer shell R < r < i?W) i-^., Ps = N^/fl^,, implying an associated charge 
density ^out = -QV^s- 

On this basis, we are able to explicitly calculate the variational free energy in virtue of Eqs. 
p2|) to ()l(i|l . As far as the entropic contributions are concerned, an analytical computation 
is feasible in a straightforward manner. For reasons of clarity, we leave out any intermediate 
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steps and present our final results only: 



N, 
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log 



LJ 



- 1, 



(18) 
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log 



(19) 
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Now, we want to investigate the electrostatic term Ues in more detail. To begin with, we 
have to derive expressions for the potential functions $in(i") and $out(r)- This computations 
are rather technical and we will outline the course of action only. In general, the electrostatic 
potential $(r) due to a charge density q{y) in a dielectric medium is given by 



e ./ r' - r 



(21) 



In order to calculate the integral above for a chopped sphere, yielding the function $in(r), 
we decompose it in infinitesimally thin discs which are oriented perpendicular to the 2;-axis 
and cover the whole sphere. The radii of these discs obviously depend on their position with 
respect to the sphere's center. Afterwards, each disc is object to further decomposition into 
concentric rings. Provided that the charge density is spherical symmetric, as is in the case 
under consideration, the charge carried by each ring can be easily calculated. The electro- 
static potential of a charged ring is known from literature , ''^1'''^ therefore the potential of a 
disc can be derived by integrating over all corresponding rings. An analytical determination 
of this integral is possible. The potential of the chopped sphere itself can then be obtained 
by another integration over all discs. The latter cannot be performed analytically, thus 
one has to resort to a simple, one-dimensional numerical integration. For a complete and 
explicit description of the procedure, we refer to Ref. 0, Appendix A, where an almost iden- 
tical derivation was put forward. So far, the potential $out (r) caused by the hollow chopped 
sphere of volume f2out containing the free counterions remains to be acquired. Thereto, we 
employ the superposition principle. The hollow region of uniform charge density fout(i") 
can be apprehended as a superposition of two chopped spheres with radii R and R^, and 
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charge densities — f?out(r) and f)out(r), respectively. In doing so, the problem is reduced to the 
calculation of the electrostatic potential of a chopped sphere with uniform charge density, 
which can be computed following the method described above for the inner sphere of volume 
fiin (see also Ref. 3, Appendix B). Clearly, with knowledge of both $in(r) and $out(r); the 
dimensionless functions h^-p and thus the electrostatic energy Ues can be obtained according 
to Eqs. (fT^ and (fT3j) using numerical techniques. 

In principle, the electrostatic-entropic interaction potential Vee{D) is obtained by adding 
up the entropic and electrostatic contributions following Eq. ()12|) and minimizing the 
Helmholtz free energy J^ee{D; R, {Ni}) with respect to R and {Ni}, cf. Eq. (fTn|) . How- 
ever, the star radii are in good approximation unaffected by the center-to- wall separation 
D, as confirmed by our simulations. The physical reason lies in the already almost complete 
stretching of the chains due to their charging. In the simulation the radius R of the star 
was measured according to 



where r^^ stands for the position vector of the last monomer of the u-th arm of the star 
and fq the core position. Fig. |21 illustrates the weak R vs. D-dependence for an exemplarily 
chosen parameter combination. Therefore, it is convenient not to determine R through the 
variational calculation, but to use average values (R)^ as obtained from MD simulations 
instead. The latter are comparable to the corresponding radii for isolated PE-stars according 
to Ref.y, see Table HI 

The amount of condensed counterions A^^i was measured by counting the number of such 
particles separated from the monomer ions by a distance smaller than Ab and performing a 
statistical average. The total number of captured counterions, A^in = N1 + N2, was measured 
by counting all counterions within a sphere having the instantaneous, arm-averaged center- 
to-end distance and again taking a time average. Since the {N^} are related through the 
equation Nc = Ni + N2 + N3, only two independent variational parameters remain, say A^i 
and N2. In our simulations, we have found that the number of condensed counterions, A'^i, 
is also approximately constant with respect to D, see Fig. El Hence, we will treat Ni as a fit 
parameter, held constant for all D and chosen in such a way to achieve optimal agreement 
with simulation results. Therefore, J-'ee{D; R, {Ni}) is minimized with respect to the number 
of trapped counterions, N2, only, reflecting the possibility of redistribution of counterions 




(22) 
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between inside/outside the stars as the distance to the wall is varied. 

Fig. lUshows a comparison of the entropic and the electrostatic contributions S and Ucs to 
the electrostatic-entropic effective potential, V^g. As one can see from the exemplary plot, the 
total entropic term 5* is the major contribution to Kc and therefore determines the functional 
form of the latter, while the electrostatic term is of minor importance. At first glance, this 
may seem to be counterintuitive in a PE-system, but electrostatics do in fact indirectly 
affect the interaction potential: Due to the presence of charges, the conformations and thus 
the radii of PE-stars are strongly changed compared to the neutral star case, leading to an 
increased range of the interaction. In the inset of Fig. ID it can be seen that total entropy of 
the system at hand is mainly determined by the trapped counterions' contribution 5*2. As 
one can see, 6*3 is only weakly influenced by the star-to- wall separation D, and 5*1 is even 
completely independent of D, as was evident from Eq. (fTHjl . We have to emphasize, however, 
that even though Si contributes a constant value only and therefore does not influence the 
effective force Fee{D) at all, the number of condensed counterions Ni nevertheless plays an 
fundamental role in our problem. Since N2 = A^in — A^i set the overall scale for of the term 
5*2, A^i becomes relevant in renormalizing the effective interaction. The treatment of A^i and 
the net charge, consequently, as a fit parameter is a common approach for charged systems, 
known as charge renormalization.i^^*^^^ 

B. The chain compression contribution 

In anticipation of Sec. IIVI we want to point out that significant deviations arise between 
the effective star-wall forces as obtained by our computer simulations and their theoreti- 
cal counterparts F^e calculated using the formalism presented in Sec. 1111 Al for intermediate 
center-to-surface distances D/R ^ 0.6 . . . 1.3. Now, we want to elucidate the physical mech- 
anism leading to these deviations and derive simple expressions for additional contributions 
to both the total effective forces and potentials, and V^. Therewith, we will conclude the 
development of our theoretical approach. 

The need to introduce these supplementary contributions is due to so far unconsidered 
conformational changes enforced by the presence of the confining wall. During the construc- 
tion of the mean-field part of our theory, we neglected such changes by assuming density 
profiles of monomers and counterions which are undisturbed compared to the case of isolated 
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stars in full 3d geometries. The influence of the wall became manifest in a truncation of the 
spheres representing star and associated counterion cloud, only. But in reality, for D in the 
order of the typical length of an arm of the star, the star will undergo strong conflgurational 
variations to avoid the wall since the monomers are not able to interpenetrate it. Thereby, 
due to the presence of neighboring arms, it can be energetically favorable for chains directing 
towards the wall to compress instead of bending away from the surface, although this com- 
pression leads to an extra cost in electrostatic energy. The latter obviously originates in a 
decrease of the ion-ion distances along the backbones of the chains. Figs.O^a) and|31^b) depict 
the situation. Once a critical value of the chain length is reached, a further shortening of the 
chains under consideration becomes disadvantageous. Electrostatic repulsions and excluded 
volume interaction increase more and more strongly, the chains will preferably curve and 
start to relax in length, see Fig. Efc). The occurrence of such compression-decompression 
processes is proved by our simulation runs. 

In what follows, we will use a simplifled picture and model the affected arms as rigid, 
uniformly charged rods of common length L and diameter (Tlj for both the compressed and 
the bent regime. In the former case, we assume the rods to have orientations perpendicular 
to the wall surface. In the latter case, we will account for the bending only via the re- 
lengthening of the rods and their change in orientation with respect to the 2;-axis. Fig. Efc) 
visualizes this approximation, whereby the shaded straight chain represents the imaginary 
rod that replaces the bent arm within the framework of our modeling. The total charge 
of one such rod is just qeaN, leading to an L-dependent linear charge density, r], given by 
1] = qeaN/L. We consider then the self-energy of one rod, as a function of its respective 
length L, to be made up of a purely repulsive contribution which is electrostatic in nature 
and arises due to like-charge repulsions, and a second harmonic term with an effective spring 
constant k^s describing the binding of the chain monomers in a coarse-grained fashion. Thus, 
we have}^ 



From Eq. (j2Sl), we then obtain the corresponding force Fi.od(-^) = —dU^od/dL that acts on 



pU,od{L)=U,,^{L)+UML) 




(23) 
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the ends of the rod parallel to its direction (a negative/positive sign of F,.od{L) corresponds 
to a force that tends to compress/stretch the rod.) The competition of a repulsive and an 
attractive part results in a finite equilibrium length Lq of the rod, i.e., Ffod vanishes for 
L = Lq. Clearly, Lq is related to the above introduced spring constant k^s by the following 
condition: 

(24) 



, 1 dUrcp 

Lo dL 

Now, let ii'min and -Dmax determine the range of center-to-surface distances for which the 
above described processes take place. Moreover, Dq denotes the critical length for which the 
transition from the compressing to the bending regime is observed. According to our MD 
results, we will fix D^^JR = 1.3, Dq/R = 0.9, and D^^JR = {2Dq - D^^^)/R = 0.5. In 
addition, we require Lq = -Dmax in what follows. We know from our simulation runs that 
it is a reasonable first-order approximation to assume the length of the affected chains as a 
function of D to be: 

/■ 

-Dmax -D G [0, -D]-Qin[ 

2Do-D De[D^,in,Do[ 



L{D) 



(25) 

D DE [Do, Dmax] 



-Dmax -D e [-Dmax, Oo[. 

By using this empirical fact, we implicitly include effects due to chain bending and entropic 
repulsions of neighboring chains, even if we did not consider corresponding energy contri- 
butions explicitly in Eq. (|23|). Therewith, a promising estimate for the chain compression 
contribution to the total effective force is: 

I otherwise. 

Here, fes is the total number of affected chains. Assuming that the chains are regularly 
attached to the colloidal core, we expect a linear relation between /eg and /, namely f^s = 
f / fo- Simulation data indicate /o = 4 to be a good choice for all parameter combinations 
under investigation. Note that the pre-factor D/L{D) = cos 7 results from geometrical 
considerations, as can be seen from Fig. Efc). 

Based on Eq. ()26|). we obtain the corresponding energy term by a simple integration: 

K{D)= r F,{D')dD'. (27) 
J 00 
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Finally, the total effective forces and interaction potentials are obtained as the sum of the 
electrostatic-entropic and compression terms: F^s = F^e + Fc and Vcs = Ke + K, respectively. 
Fig. ini shows a concluding comparison of the full effective potential V"cfr(-D) and the same 
without the inclusion of the compression term Vc{D). According to the plot, the theoretically 
predicted total potential is purely repulsive and ultra-soft. As we will see in Sec. IIVI the 
total effective forces show a striking step-like shape for intermediate distances. This is in 
contrast to the well-known case of star-star interactions.^^ 

IV. COMPARISON AND DISCUSSION 

In this section we test our theoretical model against corresponding simulation results 
to confirm its validity. Using standard MD techniques, a straightforward measurement of 
(effective) interaction potentials is not possible. Thereto, one would have to apply more 
sophisticated methods.^^^ But since the mean force acting on the center of a PE-star can 
be easily received from computer experiments,—"^ we will focus on a comparison of such 
effective forces and corresponding theoretical predictions. To be more precise: when choosing 
the colloidal center of a star at a fixed position rcore as effective coordinate in our simulations, 
the effective force can be measured as the time average over all instantaneous forces fcore 
acting on the core (cf. also the simulation model in Sec. IHl by means of the expression: 



In all considerations to come, we will deal with absolute values of the forces only. It is 
obvious from the chosen geometry that the effective forces are in average directed perpen- 
dicular to the confining wall, i.e., parallel to the 2;-axis. Thus, the mean values of the x- and 
y-components vanish, leading to the relation Fcs = \Fcs\ = Fes,z- For symmetry reasons, the 
effective forces and potentials depend on the 2;-component of rcore alone, whereas rcore,2 = F). 
Consequently, the forces are connected to the potential by the following simple equation: 



Therewith, predictions for the effective forces can be computed starting from theoretical re- 
sults for the corresponding potential (see Sec.lIIH Fig.Ej). This allows the desired comparison 
to MD data. 




(28) 




(29) 
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Fig. [7| shows typical simulation snapshots for two different values of the center-to-surface 
separation illustrating the conformational changes a star undergoes due to the presence 
of the confining wall. In part (b) of the figure, one may particularly note the PE-chains 
directed perpendicular to the wall. These arms are object to the compression-decompression 
mechanism described in detail in Sec. 1111 1^ and sketched in Fig. El Moreover, the snapshots 
qualitatively visualize our quantitative finding that the majority of counterions is captured 
within the interior of a PE-star, cf. Fig. El We want to point out that, compared to the bulk 
case,i^ the counterion behavior in this regard is not significantly changed by introducing 
(hard) walls. In this context, we also refer to Table jll where the number of condensed 
counterions A^^i as obtained from our MD runs is shown for different parameter combinations 
together with corresponding values for isolated PE-stars in unconfined geometries. 

Fig. IHl is the core piece of this work and finally presents both our simulational and the- 
oretical results for the effective star-wall forces. As can be seen from the plots, we find 
a very good agreement for all combinations of arm numbers and charging fractions. The 
theory lines exhibit a remarkable step-like shape for intermediate values of the center-to- 
surface distance D arising due to the introduction of the chain compression term F^. The 
latter had to be added to the electrostatic-entropic forces Fee to account for conformational 
changes induced by the spatial proximity of the star to the planar wall. In doing so, the 
functional form of the theoretically predicted forces accurately renders the results found by 
our computer experiments. When using the respective star's radius as basic length scale 
for the plots, the position of the step is independent of the parameters considered. In this 
sense, the effect is universal and there is no influence by specific details of the system. Note 
that the distinct kink is an artifact resulting from the special construction of Fc, it does not 
have any physical meaning. Again, it should be emphasized that the described behavior is 
in contrast to the well-known star-star case.^"^ The difference between both systems mainly 
lies in the fact that a star is strictly forbidden to access the region of the wall, while it is 
allowed to interpenetrate another star-branched macromolecule to a certain degree. To put 
it in other words: a wall is impenetrable, while a second star is a diffuse, soft object. Here, 
we also want to emphasize that the step is no hysterisis/metastability effect, as in our MD 
simulations the center-to-surface separation D is varied by pulling the stars away from the 
wall instead of pushing them towards it. This procedure prevents any coincidental trapping 
of individual PE-chains. 
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The particular values of the fit parameter A^i used to calculate the theoretical curves in 
Fig. ISl are given in the corresponding legend boxes. Obviously, A^i significantly grows with 
both the charging fraction a and the arm number /, as is physically reasonable. Table E] 
compares the fit values used to their counterparts obtained by our simulations. Thereby, one 
recognizes systematic deviations, i.e., the fit values are always slightly off from the numbers 
computed using the MD method. Both quantities are of the order of 50-60% of the total 
counterion number, for all functionalities / and charging fractions a. The discrepancies 
emerge due to the different meanings of the quantities star radius R and number of con- 
densed counterions A^^i. Within the scope of the simplified theoretical model, a PE-star is 
a spherical object of well-defined radius, while we observe permanent conformational fluc- 
tuations of the simulated stars. Thus, the MD radius does not identify a sharp boundary, 
but determines a typical length scale only. In this sense, the star can be viewed as a fluffy 
object. Most of the time, there are monomers, monomer ions and condensed counterions 
located at positions outside the imaginary sphere of radius i?, and such counterions are not 
counted when measuring A^^i, thus the theory value typically exceeds the simulation one. 
With increasing arm number / and charging fraction a, there is less variation in length of 
different arms of a star, resulting in smaller deviations between the theoretically predicted 
A^^i and its MD equivalent. Table |l] confirms this trend. 

A last remark pertains to the behavior for small center-to-surface distances D. Since our 
simulation model includes a colloidal core of finite radius i?d, the forces obtained from MD 
runs diverge for vanishing wall separations due to the core-wall LJ repulsion which mimics 
the excluded volume interaction. In contrast, our theoretical model has no core, for this 
reason the predicted forces remain finite for the whole range of distances. According to this, 
our theoretical approach is in principle not capable of reproducing a divergence of the forces, 
even in case of close proximity to the confining wall. 

V. CONCLUSIONS AND OUTLOOK 

We have measured by means of Molecular Dynamics simulations and analyzed theoret- 
ically the effective forces emerging when multiarm polyelectrolyte stars approach neutral, 
impenetrable walls. The forces have the typical range of the star radius, since osmotic PE- 
stars reabsorb the vast majority of the counterions and are thereby almost electroneutral; 
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longer-range forces that could arise due to the deformation of the diffuse counterion layer 
outside the corona radius are very weak, due to the small population of the free counterion 
species. The dominant mechanism that gives rise to the soft repulsion is the entropy of the 
absorbed counterions and the reduction of the space available to them due to the presence 
of the impenetrable wall. At the same time, a novel, additional mechanism is at work, 
which stems from the compression of a fraction of the star chains against the hard wall. For 
star-wall separations that are not too different from the equilibrium star radius, the chain 
compression gives rise to an additional repulsive contribution to the force. At close star-wall 
approaches, the compressed chains 'slip away' to the side, orient themselves parallel to or 
away from the wall and thus the compression process ceases and the additional contribution 
to the force vanishes. 

The compression force could play a decisive role in influencing the cross-interaction be- 
tween PE-stars and spherical, hard colloids of larger diameter. Indeed, for this case, the 
cross-interaction could be calculated by using the results of the present work as a starting 
point and performing a Derjaguin approximation. It will be interesting to see what effects 
the cross-interaction has in the structural and phase behavior of such PE-star-coUoidal mix- 
tures. This issue, along with the study of the effects of charged walls and wall-particle 
dispersion forces, are the subject of ongoing work. 
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TABLES 



TABLE L Conformational properties as obtained from MD simulations and corresponding fit 
parameters used in our theoretical approach for different arm numbers and charging fractions. In 
addition, results from Ref. ^ are presented for comparison. The chain length is fixed to = 50. 
For our data, the Wigner-Seitz radius is -Rw/^lj = 74.44. Note that there are in part insignificant 
discrepancies between the parameters used here and in Ref. 
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FIGURE CAPTIONS 



FIG. 1: Sketch visualizing the physical situation at hand, namely showing a PE-star of 
typical spatial extent R, and its counterion background of radius Rw- The three possible 
counterion states are illustrated: Condensed, trapped, and free. For further explanations, 
see main text. 

FIG. 2: Star radii vs. the star-wall separation D as obtained by MD simulations, referring 
to parameters f = 10, a = 1/3. The mean value, averaged with respect to D, is shown as a 
horizontal line for comparison. 

FIG. 3: Fractions of condensed counterions Ni/Nf. and captured counterions N-^^/Nc — 
{Ni + N2)/Nc for PE-stars with / = 10 arms and various values of a. As one can see from 
the plots, both quantities are in very good approximation D-independent and the majority 
of counterions is captured within the interior of the star. 

FIG. 4: Comparison of the entropic and the electrostatic contributions S and Ucs to 
the electrostatic-entropic part of the effective potential, \4e, exemplarily shown for / = 10, 
a — 1/3, and A^i = 107. Obviously, the entropic contribution dominantly determines the 
/^-dependence of V^e- In the inset, S is decomposed into the different counterion species' 
contributions Si, illustrating that the term 5*2 governs the functional form of the total 
entropy. Note that the different contributions in both the main plot and the inset were 
shifted by constant values to enhance the readability of the plot. 

FIG. 5: Schematic illustrating the physical mechanism leading to the necessity to intro- 
duce the energy contribution V^. Hollow circles denote neutral monomers, gray balls are 
monomer ions. The bigger black ball is the colloidal core of the star. Counterions are omit- 
ted in the depiction for reasons of clarity. In part (c) of the figure, the dashed rod replacing 
the chain bent to the left is shown mirror-reflect cd to the right to avoid crowding. For a 
more detailed discussion of the compression effect, sec main text. 

FIG. 6: Comparison of the total and the electrostatic-entropic contributions to the effec- 
tive potential, Ves{D) and Vee{D), for / = 10, a = 1/3, and A^i = 107. The inset shows the 
functional form of the compression term Vc{D). The main text provides an explanation for 
the need of the latter contribution and contains an explicit derivation. 

FIG. 7: Simulation snapshots of a PE-star with functionality / = 18 and charging fraction 
a = 1/3. The center-to-wall distances are D/R = 0.107 (a) and D/R = 1.07 (b). Bright 
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gray balls are neutral monomers, while the dark spheres along the chains indicate monomer 
ions. The counterions are the small, dark gray spheres surrounding the star. 

FIG. 8: Effective star- wall forces with the colloidal core of the respective stars taken as 
effective coordinate, obtained from both computer simulations and our theoretical approach. 
The fit parameters A^i used are specified in the legend boxes. Here, we show results for 
stars with functionalities / = 10 (first row), / = 18 (second row), and / = 30 (third 
row). In all cases, we are considering charging fractions a — 1/h (first column), a = 1/4 
(second column), and a = 1/3 (third column). In the rightmost picture of the first row, 
we additionally included a theoretic line calculated without taking the compression term 
into account to illustrate the need of such a contribution. Since the theoretical model has, 
in contrast to the simulation model, no (hard) colloidal core, all simulation data have to be 
displaced by the core radius R^- 



26 



iVi.n, 




FIG. 1: Konieczny, Likos 



27 



1 < 1 < 1 < 

5 

if 5 i 


o R 

_<R>„ 

i 




1 . 1 . 1 . 1 



0.5 1 1.5 2 



FIG. 2: Konieczny, Likos 



28 



0.8- 



o a=l/5 
□ a=l/4 
» a=l/3 



0.6- 



0.4 



o oOoooooooo o o 



0.2 



0.8 



0.6 



0.4 



0.2 



°°anannn □ □ □ 
oooooooo o o 



0.5 



1.5 



2.5 



0.5 



1.5 



o a=l/5 
n a=l/4 
» a=l/3 



2.5 



FIG. 3: Konieczny, Likos 



29 




30 




31 




FIG. 6: Konieczny, Likos 



32 




FIG. 7: Konieczny, Likos 
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